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Non-perturbative preheating decay of post-inflationary condensates often results in a high density, 
low momenta, non-thermal gas. In the case where the non-perturbative classical evolution also leads 
to Q-balls, this effect shields them from instant dissociation, and may radically change the thermal 
history of the universe. For example, in a large class of inflationary scenarios, motivated by the 
MSSM and its embedding in string theory, the reheat temperature changes by a multiplicative factor 
of 10^2. 



Introduction. One of the dominant paradigms of 
modern cosmology is that of inflation. In recent years 
much emphasis has been given to the epoch immediately 
after inflation, in an effort to connect it to the Stan- 
dard Model (SM) or its possible extensions, such as the 
MSSM. Although this goal has not been achieved yet, 
much progress has been made in understanding the rele- 
vant processes of reheating and preheating. 

While instabilities of the condensate may lead to in- 
tensive production of particles P, 01 j it often also gener- 
ates non-topological solitons known as Q-balls (21 . Such 
solitons are prevalent in most extensions of the SM, and 
in particular in supersymmetric ones. For example, as 
noted by (5i|,||6|, the existence of Q-ball solutions is generic 
in the MSSM parameter space. Furthermore, due to in- 
flation, we expect that some of the condensates of these 
scalar fields will be displaced far from the global mini- 
mum of their potential, thus generating the right initial 
conditions for the formation of Q-balls. The relevance 
of the latter to cosmology depends on how long these 
Q-balls live. 

Much like the approach to thermalization in preheat- 
ing, Q-balls are generated by dynamical instabilities act- 
ing on the condensate. In fact, the point of view that 
we advocate in the current paper is that Q-ball forma- 
tion should be thought of as a more controllable example 
of preheating in that the Lyapunov exponents around 
the coherent condensate configuration could be smaller. 
A related fact is that part of the systems falls into a 
new stable configuration which is close by in phase space, 
namely Q-balls. As we point out, yet another similarity 
is revealed when not all the energy of the condensate ends 
up in the Q-balls. This is the generic case which is most 
relevant to cosmology. In this case scaling arguments, 
supported by simulations of Q-ball formation [ll|, indi- 
cate that the remaining condensate go into a cold gas of 
low momenta particles with large occupation number 0. 

This cold gas can decay into hot plasma, which in turn 
can bombard and destroy the Q-balls. In the current pa- 
per we discuss the details of this process and how the spe- 
cial off-equilibrium features of the gas, such as its large 



number density, can actually stabilize Q-balls at a large 
range of charges (for a review of the standard computa- 
tions of the Q-ball formation/destruction see |^). 

The basic observation is that the finite density induces 
large masses to all fields to which it couples, which are 
much larger than termal masses for the same energy den- 
sity. Thus, the gas's life-time is significantly increased, 
thereby avoiding dissociation of the Q-balls. Further- 
more, this finite density-induced large mass can also sup- 
press the direct decay rate of Q-balls. The upshot of 
all this is that the longevity of Q-balls is generically in- 
creased by the finite density effects. It should be stressed 
that finite density effects, although similar in character, 
are not the same as the effects of non- vanishing vacuum 
expectation value. 

Relying on this observation and its consequence of a 
low reheat temperature, one can for example rule out flat 
directions in modular inflationary scenarios with gravity 
mediated SUSY breaking, in which the flat directions 
have a displacement of order Mpi from the minimum. 
The use of such moduli is very common, and is almost 
inevitable for example, in models of hybrid inflation 
with gravity mediated SUSY breaking. As detailed in 
this leads to a unique MSSM flat direction, HuHd, which 
is not carrying B — L charge and can be involved in mod- 
ular inflation. This information is extremely useful for 
building any model of inflation that involves the MSSM 
flelds or those of other supersymmetric extensions. 

The purpose of current letter is to consider the effects 
of the large occupation number of the gas on Q-balls. In 
addition to the case in , in which Q-giants with charge 
Q ~ 10^* were considered, we show that the longevity 
of the Q-balls is modified all the way down to charges of 
10^^, where the effect becomes of order 1. 

High density eflFects and MSSM inflation. The 
scenario that we will focus on is the case of thick wall Q- 
balls forming from a fragmenting condensate of a scalar 
field (f>. For example, this will be the case in gravity me- 
diated SUSY breaking with the U(l) invariant potential 
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Q-balls form for K < 0, and the Q-ball solution is given 



— (po exp[— r^/i?^] expliivt], 



(2) 



where its radius, frequency, mass, and U{1) charge are 
given by 
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Hence, when a condensate with a large VEV, </) ^ to^, 
breaks up into Q-balls with a core VEV of order 
Q-balls have a large charge. 

In the case that most of the energy in the universe 
is stored in the fragmenting condensate, as in the post- 
inflationary scenario of 0, a significant fraction of the 
charge and energy ends up inside the Q-balls 0|. The 
energy per unit charge is determined by the eccentricity 
of the orbit. It is if the eccentricity is zero, and larger 
otherwise which is the generic case. If the eccentricity is 
smaller than one, but of order one, all the charge falls into 
Q-balls (which can subsequently be eroded by coUisions), 
and the remaining energy will form a cold, high density 
gas. The number density in the gas, denoted by 5(j), is 
approximately ni ~ m^d^c/jp). Note that this density 
can be quite large since (l^*/*^!) ~ (j^Q despite the fact the 
{/(I) charge outside of the Q-balls is negligible compared 
to the charge inside. 

This scenario is verified by simulations 0|, but can 
be understood by the following scaling argument. Due 
to the attractive forces in the condensate, the inhomo- 
geneities grow like^ exp(y|X|Pi), for < ^2\K\m^. 
Q-balls form when the shortest wave length become non- 
linear. At this point, this is the smallest wave-length 
scale in the non-linear problem and hence it also deter- 
mines the shortest wavelength of the gas which results 
from the same non-linear dynamics. 

This result seems to be ubiquitous to non-perturbative 
preheating processes Both in parametric resonance 

and in tachyonic preheating, the non-linear processes are 
characterized by a low UV cut-off on the modes that are 
excited. In parametric resonance, the cut-off originates 
from the frequency of oscillation of the driving field, and 
in the tachyonic preheating it is set by the negative mass- 
squared of the tachyon. In the next section we will discuss 
the approach to equilibrium in our model. 

Had there been hot plasma between the Q-balls, it 
would have bombarded the Q-balls, leading to their rapid 



decay by dissociation. For a thermal gas between the 
Q-balls this would have been the case. For example a 
thermal gas originating at 0o ^ O.lMpi (and ^ 10'^ 
GeV) would have a temperature of 10^ GeV which would 
dissociate the Q-balls rapidly. 

This is where the high density effects become impor- 
tant, as they prevent the cold gas outside the Q-balls 
from decaying instantly due to its large density. Indeed, 
this large density contributes to the masses of all particles 
that 4> couples to. Collectively denoting these particles 
by V, their effective mass is given by. 
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Note that this induced mass is due to the large variation 
of the field, despite the fact that (50) = 0. Thus the 
decay of the gas is mediated by fields of large mass and 
can therefore be estimated by 
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(again although ((50) — 0). This leads to a much colder 
plasma outside the Q-balls than the naive estimate. In- 
deed, one finds using ||SJ), that for Q-balls with charge 
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the intermediate temperature of the plasma - after the 
gas decays and before the Q-balls' final decay and subse- 
quent reheat of the universe - is given by 
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Smaller Q-balls are not generated since the decay rates 
of both the gas and of the Q-balls themselves are smaller 
than the expansion rate, rendering such objects unstable. 

On the other hand, the temperature above which dis- 
sociation takes place was found to be 0, 0| , 
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Hence, contrary to naive expectations, for Q-balls with 
charge ^ 10"'^'', amounting to the majority of stable Q- 
balls in the spectrum, the decay of the gas cannot cause 
dissociation of these objects. We stress, once more, that 
this result is due to the high density of the non-thermal 
gas. We also note that in the results presented here we 
have neglected collisions of Q-balls, which do not change 
the above conclusions. For the case of Q-giants this was 
shown in 0, but it is true more generally. 

In particular, Q-giants, with charge of order 10^* are 
not dissociated and survive well into nucleosynthesis. As 
shown below, the reheat temperature after the final decay 
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of the Q-giants is given by0 
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In this case, the high density effects results in a reheat 
temperature twelve orders of magnitude smaller! 

As promised in the introduction, in the case of Q- 
giants generated along the flat directions of the MSSM, 
such considerations exclude all flat directions with K <Q 
from having displacement of order Mp. One can further 
exclude K > Q flat directions with lepton or baryon num- 
ber. This leads to a unique viable flat direction in the 
MSSM, H^Hfi, pinpointing the inflationary sector in the 
MSSM. 

The approach to thermal equiUbrium. Once 
formed, the non-thermal gas slowly approaches thermal 
equilibrium. In some respects, the relevant processes are 
similar to those discussed in thermalization of the usual 
preheating scenario 12j, however, there are some crucial 
differences. One such difference is that during the forma- 
tion of Q-balls, no large number densities of fields other 
than (j) are present. In such fields are initially pop- 
ulated via rapid parametric resonance or tachyonic pre- 
heating. In our case, the VEV gives mass to a heavy 
particle X (which in turn may couple to some particles) . 
For generic trajectories of the complex VEV, with ec- 
centricity of order one, the mass of X is much larger than 
the period of oscillation of (j) (of order m^) and hence 
cannot be excited by parametric resonance. Also, in 
there is no field that can undergo tachyonic preheating. 

At large occupation number, the number density of a 
thermalized gas behaves as rik oc 1/k for k ^ m and 
Uk — const for k <^ m. In our case, initially the gas is 
highly populated at momenta fc, ~ \/Km^ <C and so 
it's spectrum is sharply cut off at fc,. During thermaliza- 
tion two distinct processes occur: (i) At low momenta, 
the number density evolve to become constant as in the 
thermal distribution, (ii) Higher momentum modes be- 
come populated. As in (l), [l^, the rate of the former is 
enhanced due to the large occupation number at low mo- 
menta. Thus the relaxation time for low momentum is 
small. The production of high momentum modes, how- 
ever, must involve higher order mode-mode interactions 
(i.e. strong turbulence), which are not enhanced since the 
high momentum modes have low occupation number. 

Relying on the these observations, we can integrate 
out X and obtain that the relevant decay rate for these 
many-body processes is suppressed by the mediation of 
heavy particles and can be estimated by, 
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This estimate is not limited to a 2 — > 2 (which actu- 



ally hardly contributes to populating high momentum 
modes), but is true for all n-body interactions for the 
specific potential Hence we expect it to be valid in 
the large occupation number case at hand. 

Eq. ifT^ is of the order the decay rate of the gas, 
eq. (|SJ|. Hence, while thermalization of low momentum 
modes occurs quickly, the gas never reaches true thermal 
equilibrium and decays before it has time to populate 
high momentum modes. Had thermalization occurred 
faster, considerably higher temperature and consequent 
dissociation of the Q-balls would have resulted. 

Q-ball decay. The computation of the decay rate 
of the Q-balls also require modification in the presence 
of the non-thermal gas (although the effect is much less 
pronounced than the elimination of the high temperature 
dissociation process as discussed above). It also clarifies 
some of the physics involved. Hence below we reproduce 
the corrected decay rate and reheat temperature of the 
Q-balls, taking into account the effects of the gas. We 
will mostly follow the techniques of Ref . [l3| . 

In the core of the Q-ball the VEV of (j) is large, the 
masses of X are large and the decay is suppressed. This 
is the familiar computation |l0|, lldj llJI ■ However, outside 
of the core, the novel feature is that due to the high den- 
sity of the gas, X is massive even in the region between 
the Q-balls where the coherent VEV {(j)) vanishes, again 
suppressing the decay through the surface of the Q-ball. 
In the usual computation, one defines the core radius Tc 
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such that the decay rates are 

r(r <r,) ^m3/8^|0(r)|2, 
r(r > Tc) g^m^/STT. 



(15) 
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In our case, there is a different time dependent scale 
rd which replaces Vc- The reason for this is that in a 
mixed Q-ball/cold gas system the mass of X is m\ ~ 
5^(1 + The latter contribution is given by 

(j)Qf{t){ai/a)^, where f{t) takes into account the decay of 
the gas into light particles and mode-mode interactions, 
both of which reduce (|(5(/)p)gas. Following the discussion 
above, it may be approximated by 1 until the gas decays 
after which it rapidly approaches 0. We define rd to be 
the location where these two contributions are equal 



Td = -Rmin 



■ln((a/a03//5'),7c =7rfi?. (17) 



The decay rates are now given by 

r(r<r,(<)) ^m3/87r|0(r)p, (18) 

r(r > rd{t)) ^ g^ml/8niml + ^lf{t){ajaf).il9) 

The above has the correct limit, / ^ 0, in the case the 
gas quickly decays and one returns to the usual vacuum 
decay rate of the Q-balls. 
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The reheating temperature at the end of Q-ball decay 
can be evaluated as follows. The effective decay rate of 
the Q-ball is obtained from the charge depletion rate, 



Tn = - 
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where we have only written the leading order term, given 
that 7rf > 1. This leading contribution comes from the 
region r < ra- To this end the suppression of the decay 
rate relative to the usual scenario occurs because 7d < 7c- 
Since however, 7^ is logarithmic in time, this suppression 
is only mild. One way to see this is by solving the result- 
ing Boltzmann equations 0: 



dp^ 
dtfg 



dQjQ{t,Q)TQ{t,Q)M{Q) (21) 
QYQ{t.Q)dQfQ{t,Q) (22) 



where is the radiation energy density from the Q-ball 
decay, M{Q) is the mass of the Q-ball with charge Q, and 
fQ{tiQ) is the Q-ball charge distribution function. This 
was done in 0, where the solution to these equations 
was presented in for case of Q-giants, without, however, 
accounting for the new decay rates. One may consistently 
treat 7^ as constant, and extract the time of decay by 
comparing Tq to the Hubble scale H. One finds 
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where ti is the time at which the Q-balls form. This 
gives only a suppression to Tq relative to the conven- 
tional computation, of order jj^/jc — 1/5 at the time 
of decay. Plugging the above 7^ in l|20|l one obtains the 
reheat temperature, eq. H12|) . 

Conclusions. In recent years, a theory of preheating 
has emerged in an effort to understand universal features 
in the transition between the phase of a coherent conden- 
sate right after inflation, and the phase of thermal SM 
particles. In this letter we explored aspects of this tran- 
sition in a case where Q-balls dominate an intermediate 
stage between the two phases. 

We outlined key features of the non-thermal, high oc- 
cupation number gas that remain between Q-ball during 
their non-linear formation process. We explained how 
such effects can change significantly the Q-ball decay dy- 
namics. Compared to the situation with a thermal gas 
between the Q-balls, this high density leads to a suppres- 
sion of up to lO"'^^ in the reheat temperature. 

This, for example, occurs in inflationary models where 
inflation takes place at large VEVs of scalars present in 



SUSY extensions of the SM, such as the MSSM. In such 
cases, the Q-balls can attain charges of order 10^*. The 
non-thermal aspects of the gas suppress dissociation pro- 
cesses, leading to a reheat temperature below the nucle- 
osynthesis temperature. Under realistic circumstances, 
this picks out a unique MSSM flat direction H^Hd to be 
part of the inflaton dynamics. 

These finite density effects change the reheat temper- 
ature for much smaller Q-balls as well, all the way down 
to a charge of 10^**. 
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